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We study the probability distribution P(A) of the cosmological constant A in a specific set of 
KKLT type models of supersymmetric IIB vacua. P(A) is divergent at A = 0" and the likely value 
of A drops exponentially as the number of complex structure moduli h 2 ' 1 increases. Also, owing 
to the hierarchical and approximate no-scale structure, the probability of having a positive Hessian 
(mass squared matrix) approaches unity as h 2 ' 1 increases. 



INTRODUCTION 

Since string theory has many solutions, we are inter- 
ested in the probability distribution P(A) of the cosmo- 
logical constant A of the solutions. In this paper, we 
consider the supersymmetric vacua of the KKLT type in 
Type IIB string theory [HE]- To be specific, the model 
of interest in this paper is given by (setting Planck scale 
M P = 1) 

V =e K (\D T W\ 2 + \D S W\ 2 + \D U% W\ 2 - Z\W\ 2 ) , 
K =3 In (T + f ) - In (S + S) - ^ In (U l + Ui) , 
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W =W + Ae 
Wo =cj + b iUi - I c 2 + I S, 



(1) 



where T is the Kahler modulus, 5 is the dilation and Ui 
are the complex structure moduli, while flux parameters 
bi , Cj and di vary as the quantized fluxes vary [3]-[S] . Va- 
lidity of the weak coupling approximation requires that 
s = Re S > 1 while reality of the Kahler potential K 
requires that Ui = Re Ui > 0. The second term (with 
parameters a > and A) in the superpotential W is 
non-perturbative. The Wq used here for the complex 
structure sector is motivated by the toroidal orbifolds of 
T 6 [i[7]. 

After finding the supersymmetric vacuum solution, we 
can express A in terms of the parameters. Treating the 
parameters as random variables with some uniform (or 
smooth) distributions, we find P (A). 

• -P(A) diverges at A = 0~ for the single Kahler mod- 
ulus case : P(A) ~ (In |A|) 2 /2e 3 / 2 ^J\A\, where only 
the parameter c\ is treated as a random variable. 

• -P(A) peaks at A = 0~ for the model |l]) with multi- 
complex structure moduli. In one scenario, the ex- 
pected value (|A|) of A is shown in FIG. [I] as a 
function of h 2 ' 1 (note that ft. 2 ' 1 ~ 0(100) in many 
known models). Because of the long tail in P(A), 
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FIG. 1: The comparison of (|A|) (red circle), |A| 80 % (blue 
square), |A| 50 % (purple diamond), |A| 10 % (green triangle) are 
shown as functions of the number of complex structure moduli 



ft 2 ' 1 = 1,5,10,15,20,25,30. 



we see that (|A|) does not drop as h 2 ' 1 increases. 
To get a better feeling of the property of -P(A), let 
us introduce |A|y%, defined by /_i a i P(A)dA = 
Y%. That is, there is a Y% probability that A will 
fall in the range < |A| < |A| Y% . FIG.[l]showsthe 
result of a specific uniform distribution of the flux 
parameters. Here, typical |A|y% decreases expo- 
nentially as h 2,1 increases. For the sake of discus- 
sion, if we take |A| 50 % to be the likely value of |A| 
and assume its behavior in FIG. |l] extends to larger 
ft 2 - 1 , we see that |A| ~ 10-0.82^+2.7 for ^2,1 > 5 ^ 

• Another interesting property is the positivity of the 
Hessian (didjV). Here we see that the probabil- 
ity of having a positive Hessian increases quickly 
towards unity as h 2,1 increases. This property is 
due to the approximate no-scale behavior and the 
hierarchical structure adopted. The latter is fully 
justified for large h 2 ' 1 . 

This preference for a vanishingly small A is interest- 
ing since energetics alone may suggest more negative As 
should be preferred. This property is a consequence of 
the probability theory for the specific functional form 
of A in terms of the flux parameters [S]. This result 
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strengthens our earlier observation for positive A, where 
the functional form of A is only known approximately 
[3 [9] . That the Hessian is mostly positive allows us to 
be optimistic in the search for de-Sitter vacua in KKLT 
type models. We shall comment on what needs to be done 
to find a peaking P(A) when supersymmetry is broken 
in KKLT models. 



THE SINGLE KAHLER MODULUS CASE 

We first focus on the simple model for a single Kahler 
modulus stabilization at supersymmetric vacua as in [T], 



K K = -3 In (T + T) , W = W + Ae 



-aT 



(2) 



where Wo — c\ here. Now we solve the supersymmetric 
condition DtW — drW + (dxK)W — 0. The imaginary 
part of T has a cosine type of potential and therefore the 
solution stays at Im T = 0. Then the condition is written 
by 



-3y = (2x + 3)e~ 



(3) 



where we defined x = at = a Re T and the parameter 
y = Wo /A. The real part of the modulus is defined to 
be positive so that the volume of compactification is cor- 
rectly measured. Therefore the RHS of |3| is not only 
positive, but has an upper bound value of 3. The corre- 
sponding condition for the LHS implies — 1 < y < 0. If 
we rewrite in terms of X = —x — 3/2, ^ becomes 

:! ; ' Xe x . (4) 
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where z ~ 0.3?/. The solution of this equation is known as 
a Lambert W function or a product logarithm. There are 
two branches of solutions : W-i for X < — 1 and Wo for 
X > — 1. Since t measures the volume of the 6-manifold 
for compactification, we need X < —3/2. This branch of 
solution can be expanded around z < by 



X(y)=W- 1 (z) = L 1 (z)-L 2 (z) + 



(5) 



where L\{z) — ln(— z) and L^{z) = ln[— ln(— z)\. So 
x = —(X + 3/2) is solved in terms of the parameter y. 
Inserting this into the potential V = e Kl< (Kj?\D T W\ 2 — 
3|V^| 2 ), we get 



A = V| min = -3e K *|Wf 
=a s A 2 e 3+2W -^/[9 + 6W-i(z)} 



(6) 



81n(-y) 

We see that small A is precisely related to small y or Wo- 
The stability of the vacuum can be understood easily. 
Using the condition ([3| to eliminate A, we obtain 



Since x > and a > 0, the stability condition is 
automatically satisfied for this model 

The Probability Distribution of A - The probability 
distribution P(A) is easy to obtain : 



P(A) = J dyP 



a S A 2 e 3+2X 

9 + 6AT 



-A 



(8) 



where P(y) is the distribution of — 1 < y < and 
X(y) = yV_i(3e" 3/2 2//2). If we set A = 1 for simplicity, y 
obeys the same distribution as that of Wq. The resulting 
probability distribution P(A) is easy to evaluate numer- 
ically. It is shown in FIG [2j where a = A = 1 and Wo 
has an uniform distribution in the range — 1 < Wq < 0. 
We see that it peaks at A = 0~. 

It is interesting to see analytically the peaking behavior 
of P(A). Given P(y), the integration in ^ can be easily 
performed using the property of delta function S(g(y)) = 
\g' {yo)] 1 ^{y ^ Vo) ■ So we need to express y as a function 
of A. That is, let us solve §6§ for y. With a = A = 1, 
let us choose a uniform distribution for — 1 < y < 1. 
Introducing p = — (3 + 2X), Q becomes 



- 1/3A =pe f 



(9) 



Now the physical constraint for X(y) = W-i(z) as in ^ 
requires p to satisfy X = — (p + 3)/2 < —3/2. So the 
solution for p > is given in terms of the other branch 
of the Lambert W function, namely, 



P =W - 
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(10) 
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In [3Aln(-3A)] 



Rewriting ^ using the basic relation of the Lambert W 
function, we have 



A = 



-e-^y 



(11) 



Since now X(y) is related to p 
after using 



X 2 (9 + 6X)' 
= W (-1/3A), we obtain, 



-W (-l/3A)/2 



Wo 



1 

3A 



(12) 



where the overall sign of y is dictated by 

Using the solution ( 12 1, we get, after some calculations, 



P(A) =i [3 + 2W-i (/(A))] 2 e -3/2-W_ 1 (/(A)) i (13) 

1 



/(A) = - 



1 



3 + Wo - 



3A 



,-3/2-W (-l/3A)/2 



Using the expansion formulae of the W functions, we find 
A^o- (3 + In |A|) 2 



P(A) 



2e 3 / 2 v/fA| 



(14) 
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FIG. 2: The probability distribution P(A) of the single Kahler 
modulus model |2|. The (blue) bars are numerical result for 
Q while the curve is for the analytical formula (13 1. Both 
agree nicely even at A = 0~, where -P(A) is divergent. 



Thus we see that P(A) is actually divergent as A — s- 0~. 
We present both the "analytical result" ( 13 1 and the nu- 
merical result ^ in FIG [2] Even though P(A) is diver- 
gent at A = — , it is properly normalized in this paper, 
i.e., fP(A)dA = 1. 



SOLVING THE MODEL 

Solutions - Now we are ready to solve the model ([I]) 
with ft, 2 ' 1 number of complex structure moduli for su- 
persymmetric vacua, i.e., DtW — D$W — D\j.W = 0. 
Picking the Im T = Im S = Im P,- = solution, we have 



-1 Si - sc 2 



h 2 ' 1 -2b t - sd, 



(h 



2.1 



E 



ci - sc 2 " bi - sdi ' 
i—i 



3W + 2xAe- x = 0, 

2(£i + sc 2 ) ]li=i (h ~ sd i) 



W\ 



(15a) 

(15b) 
(15c) 
(15d) 



Re S, Ui = Re Pi, x = a Re T = at and ci 



Note that ( 15 3,c) form a pair of equations 



where s 
ci + Ae~ 

for s and x. Solving them determines s and x as well as 
Mi, so W^|min becomes a function of the parameters of the 
model. 

To solve this efficiently, we assume that the complex 
structure and dilation sector are stabilized at higher- 
energy scales than that of the Kahler modulus t. In 
this hierarchical setup, we first stabilize Ui and s in the 
model ([T]) and then the Kahler modulus stabilization 
can be dealt with. Under c\ = C\ + Ae~ x —> C\, 
VFlmin — > VFolmin, and (15i,b,d) reduce to those studied 
in [5] based on [7]. Although we deal with approximate 



solutions here, we can solve for all U{ even in the presence 
of the non-perturbative term Ae~ x . After determining 
Wo | min and s, we solve for x and W| m i n . The validity of 
this hierarchical setup will be checked a posteriori. 

Probability Distribution of A - Combining the pre- 
vious analysis, the A of the supersymmetric vacuum is 
now determined to be 



A= -?,e K \W\ 2 = 



1 



a 3 A 2 e 3+W. 1 (z) 



2 h2A + 1 sl[u i 9 + 6W-i{z 



(16) 



It is now straightforward to calculate probability distri- 
bution P(A) numerically, given the distributions of the 
parameters of the model 0. We see that P(A) typically 
peaks sharply at A = CP. 

To be specific, consider the following scenario. To deal 
with a divergent value of A coming from we intro- 

duce a cutoff / for flux inputs, 6, = — /, — / < < / 
while keeping —1 < Cj < 1 at each h 2,1 such that 90% of 
the A values fall within the Planck scale (Mp = 1) range, 
similarly to [5]. Note that the qualitative result will not 
change even if hi is randomized, though this will take 
more computer time. The remaining 10% of the data sit- 
ting over the Plank scale are discarded to maintain the 
validity of this supergravity approximation. As before, 
we set a = A = 1. The result is shown in FIG. [TJ We see 
that the likely |A| drops exponentially as h 2 ' 1 increases, 
even though (|A|) stays more or less constant. Similar 
qualitative properties are expected for some choices of 
distributions for the parameters, but not for others. 

What happens when we turn on a supersymmetry 
breaking term? If this term is simply an additive term 
to the potential V in ([T]) and results in an additive term 
to A, then the peaking behavior of P(A) will generically 
be erased [8] . To maintain and even enhance the peaking 
behavior in P(A), the supersymmetry breaking term 
must couple to the rest of V. So a good understanding 
of the supersymmetry breaking dynamics of the KKLT 
scenario (including back-reactions) is crucial. 

Probability of positive mass matrix - Since we con- 
sider basically supersymmetric moduli stabilization for 
both the complex structure and Kahler sectors, all eigen- 
values of the mass-squared matrix are expected to satisfy 
the Breitenlohner-Freedman bound |TD1 |TT] . However, for 
the purpose of uplifting to de-Sitter vacua (where the up- 
lifting will be exponentially small) later, it is motivating 
to investigate the probability that all of the eigenvalues 
are positive at the supersymmetric vacua. Let us con- 
sider the Hessian 



Hi.] = ^^j^lextrcmal, 



(17) 



where ipj runs through x, s,Ui. Note that the extremum 
solution is inserted after taking the derivatives. Insert- 



ing the numerical solutions into the Hessian (17), we can 
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h 2 ' 1 


1 


5 


10 


15 


20 


25 


Probability 


0.897 


0.981 


0.984 


0.989 


0.990 


0.994 



TABLE I: The probability of having a positive Hessian 



{didjV) at h 2 ' 1 = 1,5 
proaching unity as h 2 ' 1 



10,15,20,25. 
increases. 



The probability is ap- 



check the probability of positively defined mass squared 
matrix since the positivity of Hessian is the necessary and 
sufficient condition for the positivity of mass squared ma- 
trix at extremal points, due to the linear transformation. 

Now we calculate the probability of having a positive 
Hessian in our model assuming the hierarchical setup. 
The result is shown in TABLE. [I] This probability ac- 
tually increases as h 2 ' 1 increases, though the value at 
h 2 ' 1 = 10 is already very close to unity. Note that this 
probability is insensitive to the cutoff for 6j, rfj. 

The reason is actually quite simple. Before intro- 
ducing the non-perturbative term, we can estimate how 
the Hessian behaves by first restricting to the complex 
structure sector. At this stage, since we have the ex- 
act no-scale structure, the potential is simplified to be 
V = e K (\D s W\ 2 + \D Ut W\ 2 ) which is positively defined 
and bounded below. Then the supersymmetric solutions 
sit at Minkowski vacua, and the eigenvalues of Hessian 
are by definition semi-positive since the potential is con- 
vex downward at supersymmetric points. Next let us 
introduce the non-perturbative term. The correction to 
the supersymmetric solutions for the complex structure 
sector is of order of 



> 



Ae 



1 



W 



-(2x/3 + l) 2 In (-Wo/A) 



> -1, 



(18) 



where we used (J3J, ^ for small —Wo/A, and x > 0. 
Note that the Kahler modulus stabilization ([3| requires 
< —Wo I A < 1. In fact, the distribution of Wo is more 
peaked toward Wo = as ft. 2 ' 1 increases, therefore smaller- 
values of Wo become more likely [5] ■ Thus we see that the 
correction due to the non-perturbative term for Kahler 
modulus stabilization gets smaller and becomes negligi- 
ble when we increase h ' . This is the reason why the 
positivity of the Hessian for large h 2,1 is satisfied at most 
of extremal points. This also means that the hierarchical 
structure between the complex structure sector and the 
Kahler sector we have employed for effective analysis is 
actually reliable, as anticipated in pQ. 

This behavior of the probability is compatible with 
that studied in |12j . since their Gaussianly suppressed 
probability result will be applicable to our model only if 
we increase the number of Kahler moduli (lighter fields) 
to a large number. That a positive mass matrix is almost 
automatic for large h 2,1 is very encouraging in the search 
of meta-stable de-Sitter vacua in the KKLT models. 

Combining with the result in [S], a clear preference 



for a vanishingly small A, either positive or negative, 
emerges. This is a consequence of the non-trivial func- 
tional dependences of A on the flux parameters, when 
probability theory is applied to these string theory sce- 
narios. A better understanding of the distributions of the 
flux parameters will lead to a better determination of the 
likely values of A. It will be very interesting to find out 
what the likely value of A will be when we include addi- 
tional interactions into the model ([!]) or consider other 
stringy models. The present work opens the door for 
further fruitful studies. 
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